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^ A method  for  predicting  the  response  of  a wing  to  non-stationary 
buffet  loads  is  presented.  The  wing  is  treated  as  a cantilever  beam  with 
known  mass  distribution.  Using  generalized  co-ordinates,  the  vibration  of 
the  wing  is  governed  by  the  second  order  mass-spring-damper  oscillator 
equation.  The  buffet  load  on  the  wing  is  expressed  as  an  integral  of  the 
sectional  force,  which  is  a function  of  the  spanwise  location  and  time.  The 
non-stationary  load  is  represented  by  the  product  of  a deterministic  time 
function  and  a statistically  stationary  random  function.  The  time  history  of 
the  applied  load  is  segmented  into  a number  of  time  intervals.  Analytical 
expressions  for  the  mean  square  response  of  the  wing  displacement  are 
derived  using  a power  spectral  density  for  the  random  part  of  the  applied 
load,  similar  to  that  used  in  the  theory  of  isotropic  turbulence.  The  effects 
of  damping,  ratio  of  the  undamped  natural  frequency  of  the  system  to  the 
half  power  frequency  of  the  power  spectral  density,  length  of  time  segment, 
and  duration  of  applied  load  on  the  response  of  the  wing  have  been  investi- 
gated for  three  examples  of  the  load  versus  time  histories.  ^ 
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RESUME 


La  pr4sente  communication  porte  but  une  methode  do  provoke 
la  reaction  d’une  oile  soumise  t\  une  charge  non-stationnulre  do  buffeting.  On 
considere  quo  1’aile  est  une  poutre  en  porte-a-faux  dont  la  repartition  de  la 
masse  est  connue.  En  faisant  usage  a des  coordonnees  generalises,  la 
vibration  de  1’aile  est  gouvemf'je  par  Pequation  de  l’osoillateur  masse-resort 
arnortie  du  deuxieme  ordre.  La  charge  du  buffeting  sur  l’ailo  est  exprimoe 
comme  une  integrate  de  la  force  par  unit4  de  surface  de  la  section  et  qui  est 
une  fonction  du  temps  et  de  la  position  suivant  I’envorgure.  La  charge 
kistationnakc  est  represented  par  le  produit  d’une  fonction  determine  du 
temps  et  une  fonction  stationnaire  statistiquement  aUiatoke.  L ’evolution 
chronologiquc  de  la  charge  appliqueo  est  fractionnee  en  un  certain  nombre 
d’intervalles  de  temps.  On  derive  des  expressions  analytiques  de  la  reponse 
efficace  du  deplacement  de  1’ailo  en  utilisant,  comme  dans  la  theorie  de  la 
turbulence  isotropique,  une  densito  spectrale  developpec  en  puissances  pour 
la  partie  al&itoke  de  la  charge  appliqueo.  On  a otudie  pour  trois  exemples 
de  la  charge  en  fonction  de  Involution  chronoiogique  les  effets  sur  la  reponse 
de  1’aile  de  I’amortissement,  du  rapport  do  la  frequence  propre  non  arnortie 
du  systemc  but  la  frequence  ^ demi-puissance  de  la  densite  spectrale  develop- 
pee  en  puissances,  de  la  longueur  de  l’intcrvallo  de  temps  et  de  la  duree 
d ’application  de  la  charge. 
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THEORETICAL  ANALYSIS  OF  THE  TRANSIENT  RESPONSE 
OF  A WING  TO  NON-STATION  ARY  BUFFET  LOADS 

1.0  INTRODUCTION 

In  recent  years,  increased  attention  has  been  given  to  the  investigation  of  buffet  characteristics 
in  order  to  meet  the  design  requirements  of  future  aircraft.  This  arises  from  the  demand  for  greater 
maneuverability  and  gust  requirements  in  the  transonic  flight  regime.  Furthermore,  the  advent  of  the 
supercritical  wing  held  promise  for  expanding  the  buffet  onset  boundary  at  transonic  speeds. 

Historically,  the  term  buffeting  was  first  introduced  in  connection  with  the  structural  failure 
that  occurred  to  the  tail  of  a Junkers  monoplane  in  1930.  This  led  to  wind  tunnel  testing  of  tail  vibra- 
tions excited  by  the  separated  flow  of  the  wing.  The  complexities  of  the  buffet  phenomenon,  which 
involves  the  solutions  of  the  vibrations  of  the  wing  surface  and  the  unsteady  aerodynamics  around  it, 
made  the  problem  not  easily  amenable  to  theoretical  analysis.  Most  of  the  earlier  investigations  in- 
volved wind  tunnel  tests  to  predict  buffet  onset  and  to  provide  some  indications  of  buffet  intensities. 

A complete  account  of  the  earlier  work  on  the  study  of  buffeting  is  given  by  Pearcy  (Ref.  1). 

The  use  of  statistical  theory  in  the  investigation  of  buffeting  was  pioneered  by  Liepman 
(Ref.  2)  who  examined  the  problem  of  the  lift  force  exerted  on  a two-dimensional  thin  airfoil  moving 
in  turbulent  air.  Later,  he  extended  the  method  to  wings  of  finite  span  (Ref.  3).  The  analysis  was 
generalized  by  Ribner  (Ref.  4)  using  a model  of  turbulence  represented  by  the  superposition  of  plane 
sinusoidal  shear  waves  of  all  orientations  and  wavelengths.  These  analyses  assumed  the  exciting  force 
to  be  statistically  stationary,  so  that  well  known  statistical  concepts  developed  in  other  branches  of 
science  could  be  applied.  Experimental  studies  involving  model  testings  in  wind  tunnels  and  flight  tests 
were  carried  out  mainly-by  NACA  using  strain  gages  developed  for  measuring  the  wing  root  bending 
moment.  These  investigations  were  reported  in  detail  by  Huston  and  Skopinski  (Refs.  5,  6),  Huston 
(Ref.  7),  and  Skopinski  and  Huston  (Ref.  8).  After  this  period  of  intense  studies  in  the  fifties,  research 
in  this  subject  was  continued  into  the  sixties  at  a slower  pace.  Aside  from  wing  buffeting,  work  on 
dynamic  response  of  launch  vehicles  to  transonic  buffeting  was  carried  out  by  Cole  (Ref.  9),  Cole, 
Robinson  and  Gambucci  (Ref.  10),  and  Coe  (Refs.  11,  12).  In  the  seventies,  the  rapid  advances  in 
transonic  aerodynamics  generated  numerous  studies  on  the  transonic  buffet  characteristics  of  wind 
tunnel  models  as  well  as  of  actual  aircraft.  Hollingsworth  and  Cohen  (Ref.  13)  reported  comparisons 
of  flight  tests  with  wind  tunnel  data  for  the  F-4  aircraft.  Mayes,  Lores  and  Barnard  (Ref.  14)  studied 
the  transonic  buffet  characteristics  of  a 60  degrees  swept  wing  fighter  aircraft.  The  flight  and  wind 
tunnel  tests  on  the  F-8D  aircraft  was  carried  out  by  Damstrom  and  Mayes  (Ref.  15),  and  the  buffet 
studies  of  the  F-8  supercritical  wing  aircraft  was  reported  by  DeAngelis  and  Monaghan  (Ref.  16). 
Detailed  investigations  of  the  Northrop  F-5A  using  an  extensively  instrumented  aircraft  to  study  the 
dynamic  buffet  pressure  distributions  on  the  wing  surfaces  and  the  responses  during  a series  of 
transonic  maneuvers  were  reported  in  References  17  and  18  by  Hwang  and  Pi.  A first  generation 
method  for  predicting  the  buffet  loads  and  responses  was  described  by  Mullans  and  Lemley  (Ref.  19), 
while  Cunningham  et  al  (Ref.  20)  used  fluctuating  pressure  data  obtained  for  a rigid  wind  tunnel  model 
as  the  forcing  function  to  calculate  the  aircraft  dynamic  response  and  presented  results  for  the  F-111A. 
In  general,  methods  employing  wing  fluctuating  pressure  data  are  found  expensive  due  to  the  work  and 
time  involved  in  model  construction.  Butler  and  Spavins  (Ref.  21)  used  a more  direct  approach  in  pre- 
dicting buffet  response  in  flight.  The  technique  was  based  on  a method  proposed  by  Jones  (Ref.  22) 
who  suggested  using  a nondimensional  representation  of  the  overall  aerodynamic  excitation  and  damp- 
ing parameters  determined  from  model  testing  in  a wind  tunnel.  Preliminary  results  for  a combat 
trainer  aircraft  showed  good  agreement  between  predicted  and  measured  response  data. 

All  the  prediction  methods  for  the  dynamic  response  of  a wing  to  buffet  loads  assume  the 
exciting  force  to  be  statistically  stationary.  This  assumption  is  valid  if  at  a certain  buffet  condition, 
the  aircraft  flight  path  is  unchanged  within  a long  analysis  time.  In  wind  tunnel  testing,  statistical 
stationarity  requires  running  time  to  be  very  long  compared  to  some  characteristic  time,  for  example, 
the  natural  period  of  the  first  bending  mode,  assuming  the  first  bending  mode  to  bo  the  dominant  one. 
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However,  in  actual  flight  during  buffet  maneuvers,  the  condition  of  stationary  may  not  be  met,  or  in 
wind  tunnel  testing,  transients  have  to  be  considered  if  the  running  time  is  very  short  or  the  model  is 
pitched  during  a test.  To  analyse  the  response  of  a wing  under  such  conditions,  nonstationary  tech- 
niques to  predict  the  dynamic  response  has  to  be  developed.  A theory  of  non-stationary  response  of  a 
linear  dynamic  system  has  been  given  by  Caughey  (Ref.  23)  and  extended  by  Barnoski  and  Maurer 
(Ref.  24),  and  Holman  and  Hart  (Ref.  25).  Using  some  of  the  concepts  developed  in  these  earlier 
studies  a method  is  given  in  this  report  to  compute  the  non-stationary  response  of  a wing  for  specified 
power  spectral  density  and  time  history  of  tire  exciting  force.  With  this  method,  it  is  possible  to 
investigate  the  continuous  excursion  of  an  aircraft  into  the  buffe,  regime.  Information  such  as  the 
maximum  rms  response  and  the  time  delay  between  the  maximum  response  and  the  applied  force  can 
lie  computed  for  known  structural  and  inertial  properties  of  tire  wing  as  a function  of  the  duration  the 
aircraft  spent  in  the  buffet  regime. 


2.0  ANALYSIS 

2.1  Differential  Equation  of  Wing  Vibration 

Figure  1 shows  a schematic  of  the  wing  which  is  treated  as  a cantilever  beam  fixed  at  y = 0. 
The  differential  equation  of  beam  bending  relates  the  time  and  space  derivatives  of  the  displacement  z 
to  the  load  distribution  on  the  beam  which  may  include  inertia  forces,  external  or  applied  components, 
as  well  as  internal  components.  The  equation  can  be  written  us: 


♦ V(y,t) 


(1) 


Let 


*ty.t)  B w,(y)z,(t) 


(2) 


where  w(  (y)  is  the  mode  shape  function.  Here  z(y,t)  is  confined  to  a single  normal  mode  appropriate 
to  vibration  in  the  first  fundamental  bonding  mode,  and  the  subscript  “l"  is  used  to  denote  this  mode. 
Normally,  Equation  (2)  is  written  as  a summation  of  all  the  possible  modes,  but  in  almost  all  wing 
buffeting  tests,  most  of  the  energy  in  the  power  spectrum  of  the  buffet  bending  moment  is  concentrated 
in  the  vicinity  of  the  natural  frequency  of  the  first  mode.  Usually,  the  first  mode  is  well  separated  from 
the  higher  modes,  and  as  a result,  the  response  due  to  the  higher  mode  is  small  and  Equation  (2) 
suffices  in  describing  tire  vibrating  system.  To  find  z(y,t)  at  any  point  on  the  wing,  it  is  necessary  to 
determine  z,(t)  first.  When  the  wing  is  subjected  to  a time  dependent  load,  z,(t)  is  governed  by  the 
usual  simple  second  order  tn  uss-spr  uig  d am  per  oscillator  equation  in  generalized  co-ordinates. 
Equation  (1)  can  now  be  written  os: 


Mz(  + CZ|  ♦ KZ|  ■ L(t)  (3) 

where  tire  dots  denote  differentiation  with  respect  to  time. 

The  generalized  mass  M and  stiffness  K ore  defined  as: 


t> 

M - /m(y)wj(y)dy 

O 


(4) 
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b /d2w,(y)Y 

K"fHU^ry  (51 


Included  in  the  term  F(y,t)  in  Equation  (1 ) jure  two  aerodynamic  forces,  namely,  a fluctuating 
lift  force  and  a damping  force  associated  with  the  fluctuating  angle  of  attack  (z/V)  induced  by  the 
relative  velocity  of  the  wing.  If  the  fluctuating  sectional  normal  force  coefficient  cn(y,t)  is  known,  tire 
generalized  force  L(t)  can  be  evaluated  as: 


b 

L(t)  ■ q/ c(y)w,(y)ctl(y,t)dy 
0 


(6) 


while  the  generalized  damping  coefficient  C is 


c “ J / c,  (y)c(y)w|(y)dy 
* 0 


(7) 


In  the  above  equation,  e.  (y)  is  the  effective  slope  of  tire  sectional  lift  curve.  Define  the  undamped 

a 

natural  frequency  con  us: 


1 


K 

M 


(S) 


and  let 


f 


w C _ 


(9) 


Using  the  above  two  equations.  Equation  (3)  can  be  written  us 


z,  + + wnJZt  - Lft)  (10) 

M 


Equation  (7)  gives  the  aerodynamic  damping  coefficient.  When  structural  damping  is  to  be 
included,  it  can  be  added  directly  to  C to  give  tire  .otul  damping  coefficient  (Refs.  22,  26).  Another 
approach  of  introducing  structural  damping  is  discussed  by  Seaton  and  Rosenbaum  (Ref.  27)  and 
Davies  and  Womom  (Ref.  28).  The  damping  is  assumed  to  be  proportional  to  the  displacement  but  in 
phase  with  the  velocity  of  a harmonically  oscillating  system.  The  method  applies  strictly  to  harmonic 
motion  and  to  the  complex  exponential  form  of  solution.  For  random  motion,  great  caution  has  to  be 
exercised  in  using  such  a form  of  damping.  Let  the  fluctuating  sectional  force  be  written  os: 


P(y.t)  - q C(y)  c0(y,t)  w,(y) 


(11) 


.4* 


and  from  Equation  (6), 


b 

L(t)  = / F(y,t)dy 

O 


(12) 


In  the  above  equation,  the  buffet  load  L(t)  on  the  wing  is  expressed  as  an  integral  of  the  sectional  force 
expressed  as  a function  of  the  spanwise  location  and  time.  If  the  sectional  force  is  assumed  to  be  inde- 
pendent of  position  along  the  span,  the  overall  load  is  obtained  by  multiplying  it  by  the  wing  span. 
However,  the  exciting  force  thus  obtained  will  be  overestimated  since  fluctuations  at  one  chord  station 
will  not  be  in  phase  with  fluctuations  at  another  station. 


Assume  the  time  average  value  of  F(y,t)  to  vanish,  that  is,  F(y,t)  = 0.  The  auto-correlation 
function  of  L(t)  is  defined  as: 


Kl(t) 


lim  1 
T-*oo  2T 


T 


/ L(t)L(t+r)dt 
-T 


(13) 


Substituting  Equations  (11)  and  (12)  into  (IS),  RjJr)  can  now  be  expressed  us: 


Ki(r) 


b b 

/ / H,.  (y,  *y2  ;r)«ly  j dy2 

o w 


(14) 


where  Uj-ly,  ;r)  is  the  space-time  correlation  function  of  tire  sectional  force,  whose  mean  is  taken  to 
be  aero.  Define  a cross-correlation  coefficient  p(y,  ,y2;r)  such  thut 


%(yj^2-r) 


(15) 


where  F3(y,,t),  FJ(y2,t)  are  the  mean  square  values  of  Fly, ,t>  and  F(y2,t)  respectively,  Once 

U|.<y  | *y3  ;f)  is  known,  then  from  Equation  (14)  the  mean  square  value  of  the  load  L(t)  can  be  obtained 

since 


t2(t)  - Hl(0) 


(16) 


For  some  special  typos  of  flow,  H,  (y,  ,y2;r)  has  the  same  shape  at  different  values  of  y and 
depends  only  on  £ « y2  - y , . From  Equations  (14)  to  (1G), 


.. b b - y 

L2(t)  - F3(y,t)/dy/  p«WS  (17) 

o -y 


If  p($)  is  equal  to  a constant  and  has  a value  unity  over  the  interval  0 < y < b,  then 
Equation  (17)  becomes 
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L2(t)  = F2(y,t)  • b2 


which  gives  the  mean  square  value  of  the  load  for  perfect  correlation.  If  p(£)  drops  rapidly,  then 

b-y 

/ p(£)d£  hardly  depends  on  y.  Since  this  has  the  dimension  of  length,  it  can  be  termed  as  the  correla- 

-y 

tion  length  X.  From  Equation  (17), 


L2(t)  = F2(y,t)X'b  (19) 

Since  the  sectional  force  is  correlated  only  for  a short  distance,  then  X is  small  and  hence  L2(t)  is 
smaller  than  that  obtained  from  Equation  (18)  for  perfect  correlation.  In  actual  cases,  the  form  of  the 
correlation  function  RF(yj  ,y2  7)  depends  on  the  wing  geometry  (swept  or  unswept),  presence  of  lead- 
ing edge  flaps  and  other  special  characteristics  in  the  design  of  the  wing.  However,  in  the  preliminary 
stage  of  design,  it  is  helpful  to  obtain  a rough  estimate  of  the  buffet  load.  In  order  to  do  so,  some 
approximate  form  of  the  correlation  function  is  needed  since  the  exact  form  will  not  be  known  for  the 
wing  until  it  has  been  constructed  and  tested  experimentally.  Although  it  is  not  within  the  scope  of 
this  report  to  discuss  this  subject  in  great  detail,  two  other  forms  of  Rp(yj  ,y2  ;r)  are  briefly  considered 

and  the  resulting  values  of  L2(t)  are  evaluated  to  give  some  idea  of  the  form  of  the  fluctuating  load  on 
the  wing. 

For  a straight  wing  with  large  aspect  ratio,  the  influence  of  tip  vortices  on  the  fluctuating 
load  is  unimportant.  The  buffet  load  will  be  predominantly  determined  by  the  fluctuating  pressure  in 
the  separated  flow  region.  In  Reference  2,  Liepman  considered  isotropic  turbulence  to  illustrate  the 
random  excitation  on  a thin  wing  moving  through  turbulent  air.  If  it  is  assumed  that  the  buffet  load 
has  a similar  form  of  power  spectrum,  the  correlation  function  can  be  expressed  as: 


Rf(*,0)  = F2(t)  9-W 


LX 

\ 2c, 


where  c is  the  mean  chord  and  F2(t)  is  constant  and  independent  of  y.  Substituting  into  Equation  (17) 
gives 

U(t)  - F^JtHl-  e'b/c)bc  (21) 

Another  form  of  the  correlation  function  that  is  of  some  interest  is  given  by 

R,,a\0)  » F2^)  e-s*^  (22) 


where  5 is  the  spatial  decay  constant  which  is  determined  empirically.  This  form  of  the  correlation 
function  is  similar  to  that  given  by  Mullan  and  Lemley  (ltof.  19)  based  on  data  obtained  from  wind 
tunnel  tests  of  a 10%  scale  model  of  the  F-4E  aircraft.  For  highly  swept  wings,  they  found  that  wing 
tip  vortices  and  snag  vortices  have  an  important  influence  on  the  fluctuating  pressure  on  the  wing.  The 
buffet  load  for  this  correlation  function,  after  substituting  Equation  (22)  into  Equation  (17)  and  per- 
forming the  integration,  is 


(23) 


— = — — 2c 2 b 1 . . 

L2(t)  = F2(t)‘— - --  -(1-  e-6bM 
5 (e  6 


2.2  Response  to  Non-Stationary  Input 

For  a given  buffet  maneuver,  the  mean  square  loading  on  the  wing  versus  time  may  be  repre- 
sented schematically  as  in  Figure  2.  The  points  A and  B correspond  to  the  onset  and  exit  of  buffeting 
respectively.  The  duration  of  time  the  aircraft  spends  in  the  buffet  regime  is  denoted  by  TB . The  load 
is  non-stationary  and  time  varying  properties,  such  as  the  mean  square  of  the  response,  can  only  be 
determined  by  instantaneous  averaging  over  an  ensemble. 

The  load  L(t)  can  be  represented  as: 


L(t)  = 7(t)*f(t)  (24) 


where  7(t)  is  a deterministic  function  of  time  and  f(t)  is  stationary  (Ref.  29).  Furthermore,  it  is 
assumed  that  f(t)  is  Gaussian  such  that  the  response  is  Gaussian  as  well.  The  non-stationary  input  can 
be  represented  by  a time  segmentation  technique  whereby  the  time  T,,  the  aircraft  experiences  buffet- 
ing is  divided  into  a number  of  time  segments.  At  any  segment,  L(t)  is  expressed  as: 


m 

L(t)  = £ 7r(t)*fr(t)  (25) 

where  the  subscript  r denotes  the  rth  time  segmont  and  the  deterministic  function  7r(t)  may  bo  con- 
sidered as  a shaping  function  written  in  the  following  form: 


7t(t)  ■ eA(t) 


(26) 


where  0t(t)  * u(t-t,  _ j ) - u(t-t, ).  Here  u(t)  is  a unit  step  function.  et  is  a constant  at  each  time  seg- 
ment and  represents  the  intensity  of  the  generalized  input  load.  Figure  3 illustrates  the  variation  of  c 
with  time  and  nonce  represents  the  load  intensity  us  u function  of  time.  The  formulation  ns  described 
in  the  above  mariner  gives  on  approximate  description  of  a continuously  varying  load.  However,  such  a 
representation  cun  bo  exact  in  practice,  for  oxample,  in  wind  tunnel  testings  where  the  model  angle  of 
incidence  advances  in  steps  with  time  so  that  tire  buffot  load  can  be  represented  by  Equution  (25). 

To  determine  Zj(t),  tho  response  z1((t)  due  to  tire  rth  timo  segment  is  first  ovaluuted  since 


E[z,2(fc)]  - jSisfzjJtt)]  (27) 


where  m is  tire  timo  segment  where  t is  located,  thut  is,  tm„,  < t < t,,,  and  E denotes  tire  expected 
value.  Knowing  tire  load  L(t),  Equation  (10)  may  be  solved  for  ztf(t).  Introduce  a frequency  response 
function  H(oj)  as 


H(oj)  \ — •— 

M(wh2  - u>2  + iBfwtuJ 


(28) 


8 

h ■ 


■■- ■ t 


■.».> 


|H(co)|2 


M2[(wn2-  co2)2  + 4£2w2u;I12] 


The  covariance  Czill(t,  ,t2)  can  be  expressed  in  terms  of  the  power  spectral  density  of  the  applied  load 
sii(wi,w2)  as: 


" ~ t(w1trw2t2) 

a / / SLL(o>1,co2)H(co1)H1('((*;2)  e dcL>,dw2 


• . ■:  .< 
■ ( * . 


m: 


where  H*(w2)  is  the  conjugate  of  H(to2).  If  z,(t)  has  zero  mean,  then 

czi*i(ti ,t2)  c R*lzl(tt ,t2)  (31) 

where  R,,,^  ,t2)  is  the  auto-correlation  function.  The  mean  square  response  can  simply  be  obtained 
by  setting  tj“t2at,  that  is, 


E[Z|2(t)]  - Rzlzl(t,t) 


^ V * 

fry;; 

f Vv  * . . 

• •- 


[•  \ 


-•iv 


f 


In  Equation  (25),  let fr(t)  = 0.  Since  R,  L(t.,t2)  - E[L(t,  )L(t,)l , itmav  beexnressed  in  the 
following  form  after  using  Equation  (26)  for  L(tt ) and  L(t2):  't2 > J » may  be  expressed  m the 


■ rft  £ 'yr(tl)T8(t2)E[fr(tl)fs(t2)] 

Expressing  Elfr(ti)f,(tj)J  » R|(fs(tj,t2),  Equution  (33)  can  be  rewritten  as 

ID  ID 

Rll(Mj)  - E|7r(ti)7l(t2)Rflfl|(tl,t2) 

If  Ru(ti  h)  Is  known,  then  from  the  generalized  Wioner-Khinehino  relations  (Hof.  29), 


1 **  m “It Id  1 t j - Cdltl ) 

*»ll(wi,w2)  - ^5  [m  RLL(f|.t2)e  dt,dt2 


Combining  Equations  (34)  and  (35)  yields: 


• 8- 


- X.:<2 


m m oo  oo  -Kojjtj-wjtj) 

®llC^i *^2)  “ “ S S / / 7t(ti)7s(t2)Rf  f (tj,t2)  e dttdt2 

4jj.2  t«U*l-«  .«  * 


If  ft(t)  and  fs(t)  are  stationary,  then  Rf{fs(t1  ,t2)  = Rfr,s(ti  -t2)  = Rlrf&(r),  that  is,  the  corre- 
lation depends  only  on  r =>  tj-t2.  Since 


1WT 

Rfrfs(r)  = / Sfffs(cjie  do) 


and  substituting  into  Equation  (86)  gives 


in  in  00 

SLL(coi,W2)  = EE/  Sf  f (w)  A^w-Wj)  A*  (a)-w2)dco 

r°l  s»l _«  ‘ 1 


where 


1 00  Ku>-W|)t  j 

At(w-w1)  = — / 7t(ti)e  dt[ 

*7T  —00 


J «•  -i(w-w2)t2 

As\w-a2)  = — / Ts(t2)  e dt2 

#y7T  -~«*j 


From  Equation  (30),  tire  mean  square  value  of  z.  (t)  is  given  as: 


p « « » U(a)  J '"'UJJ  )t 

E lV(t)J  -/  / SLL(wl,w2)H(cu1)H,(w2)e  dcu,dco2 


Substituting  Equation  (38)  into  Equutiou  (41)  yields 


r ..  111  111  111  »> 

E U|(t)J  “ E I [».*(«]  - E E / S,  , (w)Ir(t,co)Is*(t,co)dw 

u r«  l *•  ' J r«  l »■  t -»  1 s 


where 


It(t,w)  « / Ajtw-Wt)  H(wi)  e dwt 


"m 


3,4, 


0.1  . 


vO.-V'*  *»>•*. 
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m 


I*  (t.co)  = / A*{co-w2)H*(t02)e  dw2 


I(t,co)  is  the  time  dependent  frequency  response  function.  Equation  (42)  gives  the  general  formulation 
for  amplitude  modulated  stationary  inputs  and  is  a function  of  the  spectral  density  St  r f r (co ) and  cross 
spectral  density  Sfrfs(£o). 

To  evaluate  I(t,co),  consider  the  time  segment  denoted  by  tt_i  < t < tj  (Figure  3).  It  can  be 
shown  that  in  Equation  (39) 


Ar(^i)  = ^ ^(WrW)  + ~ 


m, 

i -w)  j \ 


i(w-Wj)tr_j  j(cj  —OJ  J )tr 

e -e 


Substituting  this  expression  into  Equation  (43)  gives 


er  | iwtj  »■  H(coj)  iw((t-tf) 

= r~{e  I :e  dwj 

2 it  \ -“>  i(co~u)j ) 


iwtr_j  ■»  H(CUj ) iwi(t-tr_j) 

•e  / ~e  dwt 

Kw-ojj) 


* H(coj)  iwj(t-tr)  » H(cjj)  iu)j(t-tr_j) 

The  integrals  / e dw;  and  / e dec  ■ can  be  evaluated  by 

-«i(a>-W|)  -~i(w-cO|) 

contour  integration.  The  function  H(wj ) (w-wj  )-1  has  poles  at  aj  = oj,  o2  ~ Wj  + i f oju  and 
o3  = - + i £ to,,,  where  is  the  damped  natural  frequency  and  is  defined  as: 


<*>d  a 


Considering  the  first  integral,  the  integration  path  for  t > tr  is  taken  above  the  real  axis  as 
depicted  in  Figure  4a,  while  for  t < tr,  the  integration  is  below  the  real  axis.  Similar  paths  of  integra- 
tion are  taken  for  the  second  integral.  Upon  substitution  of  these  integrals  into  Equation  (46),  tire 
final  form  becomes: 


Iwl  K.,/  iw 

It(t,co)  » etH(w)  • u(t-t,_[)  e -e  (i/'(t-tr_l)  + — 0(t-tt_ t) 

\ 


lwt  iwt,  / iw 

u(t-tr)  e -e  (i>(t-tf)  + — 0(t-t,) 
\ wd  i 
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where 


^wnt  / COn  > 

= e (cosojdt  + f — sincj(it 
V wd  ) 


-fwnt 

<j>(t)  = e sin  wdt 


In  buffet  maneuvers  or  in  wind  tunnel  testing  of  a pitched  model  whose  incidence  changes 
with  time  like  a step  function,  the  loading  on  the  wing  can  be  approximated  by  Equation  (25)  with 
fr(t)  being  the  same  in  all  time  segments.  The  only  changes  in  L(t)  being  the  intensity.  Th;s  may  be  a 
fair  approximation  since  buffeting  is  primarily  due  to  the  fluctuating  pressure  forces  of  the  separated 
flow  on  the  upper  wing  surface.  The  frequency  structure  of  the  separated  flow  may  not  change  by  a 
large  extent  as  the  incidence  varies.  If  it  is  assumed  that  there  is  no  correlation  between  different 
segments,  then  from  Equation  (42), 


E[sit(t)]  = Strrr(w)|Ir(t,w)|2  dw 


where  after  some  algebra,  IIr(t,a))|2  is  given  as  follows: 


fort^  <t<tr, 


Ur(t,w)|2  = er2|H(w)|2  • jl  - 2i//(t-tt_j)  cos  to(t-t,._x) 


- 2 — 0(t-tr_j)  sin  oj(t-tr„ j ) 


+ ^(t-Vt)  + — -02(t-tr_,) 
Hi2 


fort>tr, 


|Ir(t,cj)|2  = et2IH(w)|2  • 02(t-tr„,)  + 02( t-tr) 


W2  T W2  1 

+ — r 02(t-tr_o  + — -02(t-tr) 

<V  wd2 


*•  v'  > 


- f.  ■ . . , 


2cosw(tr-tr_1)^(t-t,_i)^(t-tr) 


CO2  "1 

— 0(t~tr_  J )^(t— tj. ) I 

“d 


2^sinw(tr-tI_1)^(t-tr_1)0(t-tr) 


M\ 


2.3  Evaluation  of  Mean  Square  Response  for  a Given  Power  Spectral  Density  of  the  Input 

Equation  (51)  can  be  integrated  to  give  the  mean  square  response,  provided  the  form  of  the 
power  spectral  density  Sfrfr(w)  of  the  input  is  specified.  Cole  (Ref.  26)  gives  a form  of  the  spectral 
density  which  has  been  found  to  adequately  describe  the  buffeting  input  to  launch  vehicles.  In  the 
present  notations,  this  can  be  written  as 


Sffft(w)  » 


" vV 


Bfrfr(“)  o) 

where  wf  ia  the  half  power  frequency  and  S0  is  a constant.  In  Figure  5,  — is  plotted  versus  — 

S0  WI1 

<*>n  «n 

for  different  values  of  — . For  — = 0,  Sf  f (w)  =■  constant  = Su,  which  is  tire  power  spectral  density 
wf  U)f  T'r 

for  white  noise.  As  — increases,  the  energy  of  tire  exciting  force  near  the  natural  frequency  con 

(Of 

diminishes  and  this  will  have  air  appreciable  effect  on  the  response  of  the  system. 

Consider  the  rth  time  segment  and  using  Equation  (52)  or  (53)  for  |I.(t,(o)|2  together  with 
the  expression  for  |H(w)|2  given  in  Equation  (29),  the  mean  square  value  E[zli‘2(t)J  can  be  evaluated 
analytically  since  the  form  of  the  power  spectral  density  is  known.  Let 


2fM2wn 


E(zlt2(t)J 


Using  contour  integration  and  after  much  algebra  the  final  expression  for  S (zlr2(t)]  are  as  follows: 


where  for  tr_l  < t < tj, 


(67) 


.--:  o’c-  rtA  ;;r:  r^vj:^'’.  ■’c/i.  ^^'.;',- 


m 


k \ I; 
: r- 
I | 
i1  $ 

fe.  | 
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and  for  t>tj. 


£tzij(t)]-*e2 


Mt-t r-i)  + ^(t-t,) 


+ — - [«2(t-Vi)  + 02(t-t,)] 

i-r2 


- 2^(t-t*)  + 


Mt-VjMW  0 


w2 


20(t-tr) 

i-r2  |yw* 


0<t~tr_  j MW  | )-  0(t~tr„,  )0(tr-tr_l ) 


CJ„ 


(64) 


In  the  limit  for  — ° 0,  these  two  expressions  give  the  response  to  white  noise. 

Wf 

2.4  Computational  Procedure 

The  steps  involved  in  computing  the  mean  square  response  is  fairly  simple  and  straight- 
en , 

forward.  For  a given  value  of  — , S[zj,(t)J  can  be  evaluated  from  Equations  (56)  and  (57)  (or 

Wf  ’ 

Equations  (63)  and  (64)  for  white  noise)  once  the  e,’s  are  given  as  functions  of  time.  Referring  to 
Figure  3,  at  any  specified  time  t,  the  procedure  is  to  establish  the  time  segment  where  t lies.  If  it  falls 
between  t,n.j  and  then  £[zj2(t)]  is  obtained  as  follows: 


m-i 


fil*i(t)l  - filzfro(t)l  + £ S(zlf2(t)] 


(65) 


where  the  first  term  on  the  right  hand  side  of  the  equation  is  determined  from  Equation  (56)  while  tho 
second  terra  from  Equation  (57). 

Three  cases  are  considered  where  e varies  with  time  in  the  form  of  the  following  functions: 

Case  I:  : Sinusoidal 
Case  II:  Ramp 
Case  III:  Triangular. 


These  functions  are  shown  schematically  in  Figure  6 and  shall  be  referred  to  the  case  number  in  the 
future.  For  convenience,  the  maximum  value  of  e is  taken  to  be  unity  and  this  does  not  affect  the 
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generality  of  the  solutions.  The  number  of  time  segments  in  TB  will  be  specified  later  on  after  the 
effect  of  varying  the  segment  duration  (tf-tj.j ) has  been  investigated.  The  variables  involved  in  the 

analysis  are  — , £,  and  TB,  and  their  influence  on  the  mean  square  response  are  studied  in  some  detail, 
cop 


3.0  RESULTS  AND  DISCUSSIONS 

Since  Equations  (56)  and  (57)  are  the  principle  expressions  to  be  used  to  compute  the  mean 
square  value  of  Zj  (t)  for  specified  non-stationary  loads,  their  characteristics  are  first  investigated  in 
some  detail.  Equation  (56)  determines  the  response  or  rise  to  a step  modulated  stationary  exciting 
force,  while  Equation  (57)  governs  the  decay  when  the  applied  force  is  withdrawn.  Taking  er  to  be 
unity,  Equation  (56)  is  plotted  in  Figures  7 to  10  as  a function  of  time  t'-t'_j  for  three  values  of  the 

(jJn  o on 

damping  ratio  (£  * 0.02,  0.04  and  0.08)  and  four  values  of  — , namely, — = 0,0.5, 1 and  2.  The 

COf  COf 

superscript ' is  used  to  denote  that  time  is  non-dimensionalized  with  respect  to  the  undamped  natural 
period  Tn  « 2rr/con.  From  these  figures,  it  can  be  seen  that  for  larger  values  of  f , the  mean  square  of 
Zj(t')  reaches  an  asymptotic  value  for  t'-t|_j  -*•  °°  sooner  than  that  with  smaller  values  of  £.  Increasing 
the  value  of  con/cor  tends  to  diminish  the  asymptotic  value.  This  can  be  explained  from  Figure  5 where 
the  normalized  power  spectral  density  is  plotted  versus  co/con.  The  energy  in  the  exciting  force  near 
the  natural  frequency  con  diminishes  for  increasing  con/oof,  and  hence  it  can  be  expected  that  the  value 
of  the  mean  square  response  will  decrease  as  ion/oof  increases.  From  Equations  (49)  and  (50),  it  is  seen 
that  as  t'-t'_j  -*■  oo,V/(t'-t[_1)  -*■  0 and  0(t'-t4_j ) -*■  0.  Hence  Equation  (56)  can  be  written  as  (taking 

er  - 1): 
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S[zi2(t')]  decreases  as  con/cof  increases.  For  con/cof  » 1, 
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In  the  limit  con/cof  -*  «»,  &[zj2(t')]  -+•  0. 

In  the  computations,  no  significant  changes  in  S [Zf  2(t')]  compared  to  that  at  con/cof  = 0 are 
detected  for  values  of  con/cof  less  than  0.1.  These  curves  are  replotted  in  Figures  11  to  13  for  fixed  £ 
with  (Op/tOf  as  the  variable.  The  oscillations  observed  are  due  to  the  sine  and  cosine  terms  in  the 
expressions  for  i^(t'-t'_1 ) and  <f>( t'-t'_j),  but  they  are  rapidly  damped  out  for  increasing  t'~t'_ , . 

Equation  (57)  is  used  to  compute  the  decay  of  the  mean  square  of  zx  (t')  after  tire  exciting 
force  has  been  withdrawn  at  time  t' . Using  the  same  values  of  £ and  con/cof  as  before,  the  decay  to  a 
pulse  modulated  statioqary  exciting  force  are  shown  in  Figures  14  to  25  for  different  values  of  the 
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pulse  duration  (t^-t|_|  ■ 0.5, 1,  2,  4 and  8).  Since  the  curves  are  plotted  with  t -tr_  j as  the  absicca, 
then  t^-t|_i  is  simply  measured  from  0 along  the  axis  to  the  time  t,-tj_1  when  the  decay  curve  starts. 

For  large  values  of  the  pulse  duration.  Equations  (49)  and  (50)  give 


*(t;-t;.i)-o 


$(t'  - t|_1)-*'0 

as  t'  - t'_i  -*  °°. 

Also  from  Equations  (58),  (60)  and  (62), 
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Equation  (57)  can  be  written  as  (taking  er  * 1) 
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which  is  independent  of  t'-t,_| . In  other  words,  the  decay  is  independent  of  the  step  duration  for 
large  values  of  t' -tj._  j . From  the  figures,  it  is  also  observed  that  for  constant  f , the  decay  is  more  rapid 
as  cun/wf  decreases,  especially  for  cases  with  the  larger  values  of  t'-t^ _ t . This  can  be  shown  from 
Equation  (68). 

For  ojn/cof  « 1, 
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while  for  u>nj cof  » 1, 
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Also  of  interest  in  those  figures  are  the  oscillations  in  the  decay  curves.  For  small  damping, 
the  oscillations  are  very  clearly  shown.  The  amplitudes  increase  with  con/cOf  while  the  decay  is  less 
rapid.  This  is  even  more  significant  for  smaller  values  of  t^— t, _ ] and  will  be  shown  later  to  have  quite 
an  effect  on  non-stationary  loads. 

For  an  input  load  with  sinusoidal  variation  of  e with  time  (Case  I in  Figure  6),  £ [z,2(t')]  has 
been  computed  using  Equation  (65).  The  duration  of  the  applied  load  T^  is  forty  times  the  natural 
period,  and  it  is  divided  into  eighty  equal  time  segments  giving  a value  Atr  = t' -t' _ | of  0.5.  Figure  26 
shows  the  response  for  J = 0.02,  0.04  and  0.08  with  ton/cof  = 0.  Also  plotted  in  the  same  figure  is  the  e 
versus  time  curve  for  comparison  purposes.  Similar  to  Figures  7 to  10,  the  rise  time  is  seen  to  be 
smaller  for  the  larger  values  of  f.  Shown  in  the  figure  is  the  delay  r,  which  is  the  time  where  the  maxi- 
mum value  of  & [Z|  2(t')]  lags  behind  the  applied  force,  and  the  amplitude  function  5 defined  as 

5 = 1-  S[z,2(t')]max. 


The  effect  of  changing  the  number  of  time  segments  has  been  studied,  and  it  is  found  that  if 
At'  < 1,  the  differences  in  the  computed  results  are  not  very  significant.  Figure  27  shows  the  response 
curves  for  three  values  of  At'  and  it  is  seen  that  the  curves  with  At'  = 2 and  0.5  are  quite  close.  In  all 
subsequent  computations,  AtJ  is  taken  to  be  0.5. 

The  effect  of  o>n/oof  on  fi  [zj  2(t')]  are  shown  in  Figures  28  to  31  for  cun/cjf  = 0.1,  0.5,  1 
and  2 with  f = 0.2.  The  oscillations  in  the  response  curves  increase  with  con/cuf,  and  the  peak  to  peak 
values  reach  a maximum  at  con/cof  = 1.  These  oscillations  arise  mainly  from  the  second  term  of 
Equation  (65)  where  previous  results  (Figures  14  to  25)  show  that  for  increasing  values  of  con/cor, 
the  decay  to  a pulse  modulated  input  has  fairly  large  amplitude  oscillations  and  diminishes  rather 
slowly  with  time.  The  summation  in  Equation  (65)  for  a large  number  of  segments  help  to  accentuate 
this  effect,  thus  resulting  in  large  regular  oscillations. 

To  illustrate  the  effect  of  the  duration  of  the  applied  force,  results  with  TB  = 10  are  pre- 
sented keeping  all  other  variables  the  same  as  before.  Figure  32  shows  K[z,2(t')]  versus  t'  for 
wn/ojf  = 0,  f = 0.02,  0.04  and  0.08.  Comparison  with  Figure  26  shows  that  the  delay  r is  shorter 
while  5 increases.  Comparing  Figures  33  to  36  with  Figures  28  to  31  gives  the  effect  of  Tj}  on  the 
mean  square  response.  Decreasing  the  value  of  TB  results  not  only  in  a smaller  peak  value  of  S [zj2(t')J , 
but  also  decreases  the  peak  to  peak  value  of  the  oscillations. 

The  delay  r is  a function  of  f and  TB.  Its  variations  with  Tg  for  f = 0.02,  0.04  and  0.1  are 
shown  in  Figure  37.  The  delay  increases  with  TB  in  steps  of  half  the  natural  period.  For  f = 0.04, 
and  0.1,  computations  up  to  Tb  = 60  show  that  r does  not  change  from  the  values  1 and  2 which 
correspond  to  those  reached  at  T’b  = 10  and  18  respectively.  It  thus  appears  that  r has  reached  the 
maximum  for  these  two  values  of  f.  For  $ = 0.02,  computations  have  not  proceeded  for  sufficiently 
large  values  of  T'b  to  determine  the  maximum  r.  Figures  38  to  42  show  the  variations  of  8 with  Tg 
for  con/tof  = 0,  0.1,  0.5, 1 and  2 with  £ as  the  parameter.  The  results  indicate  that  8 increases  with 
decreasing  f for  any  fixed  wn/ojf  while  increasing  TB  results  in  smaller  5.  For  large  values  of  TB,  8 
tends  to  an  asymptotic  value  and  the  approach  to  this  value  is  more  rapid  as  f increases.  By  increasing 
con/cof,  5 is  also  increased  as  shown  in  Figures  43  to  46  where  5 is  plotted  against  TB  for  fixed  f 
with  0Jn/cof  as  the  parameter.  These  curves  again  indicate  that  for  larger  values  of  f,  5 approaches  its 
asymptotic  value  more  rapidly. 

Similar  to  Figure  27,  the  effect  of  At'  on  £[zt2(t')J  is  plotted  in  Figure  47  for  Case  II 
with  f = 0.08  and  three  values  of  Atj.  (0.5,  2 and  4)  are  considered.  The  same  conclusions  are  arrived  at 
as  before,  that  is,  for  AtJ.  < 1,  the  differences  in  the  response  curves  are  very  small  for  different  At', 
and  hence  a value  of  0.5  is  again  used  throughout  the  computations. 

The  results  for  Case  II  (T'b  = 20)  and  Case  III  (Tg  = 40)  are  shown  in  Figure  48  which 
serves  to  illustrate  the  effect  of  f at  con/tor  = 0.  Figures  49  to  52  give  the  response  curves  for  these 


• 18. 


two  cases.  The  value  of  f is  0.02  and  ton/t of  = 0.1,  0.5, 1 and  2.  The  effect  of  T|j  is  demonstrated  in 
Figures  53  to  57  where  T'b  = 5 for  Case  II  and  T'B  = 10  for  Case  III.  These  results  are  very  similar  to 
those  for  Case  I.  The  most  significant  differences  between  Case  I and  Case  III  are  the  smaller  values  of 
the  & [z(  2(t')lmax  for  the  latter  case.  This  is  to  be  expected  since  e reaches  its  maximum  and  decreases 
at  a much  faster  rate  than  in  Case  I. 


4.0  CONCLUSIONS 

A method  for  predicting  the  response  of  a wing  to  non-stationary  buffet  loads  has  been 
developed  and  applied  successfully  to  a number  of  hypothetical  examples.  For  buffet  maneuvers  in  an 
aircraft,  the  fluctuating  loads  change  continuously  and  hence  the  analysis  presented  herein  which 
models  the  load  by  a time  segmentation  technique,  is  to  be  treated  as  an  approximation.  However, 
this  approximation  has  been  found  to  be  fairly  good  since  a study  of  the  effect  of  the  duration  of  time 
segments  shows  that  if  it  is  below  a certain  value,  there  is  little  change  n the  results  when  smaller  time 
segments  are  taken.  In  wind  tunnel  buffet  tests,  it  is  possible  to  control  the  incidence  of  the  model  by 
pitching  it  in  a prescribed  manner.  The  buffet  load  can  be  made  to  vary  with  time  in  the  form  of  a 
step  function.  In  this  case,  the  method  of  representing  the  load  described  in  this  report  can  be  consid- 
ered to  be  exact. 

The  form  of  the  power  spectral  density  of  the  input  load  used  in  this  study  is  similar  to  that 
encountered  in  the  theory  of  isotropic  turbulence,  and  analytical  expressions  for  the  mean  square 
response  of  the  wing  displacement  has  been  derived.  The  effect  of  varying  the  ratio  of  the  undamped 
natural  frequency  to  the  half  power  frequency  wn/co,-  on  the  response  of  the  system  has  been  investi- 
gated and  it  is  shown  that  if  i oM/w,  is  small,  the  maximum  value  of  the  response  is  greater  than  that 
obtained  for  larger  values  of  wn/tu, . This  is  due  to  the  distribution  of  energy  in  the  input  load  sinco  a 
large  value  of  wn/wf  implies  that  little  energy  is  distributed  near  the  undamped  natural  frequency.  It 
should  be  noted  that  for  white  noise,  equal  energy  is  distributed  at  all  frequencies,  and  u>n/to,  in  this 
case  is  equal  to  zero. 

The  duration  of  the  applied  load  T|t  and  the  form  of  the  time  history  of  the  load  have  also 
Ireen  studied.  Detail  results  are  presented  for  a sinusoidal  variation  of  tire  force  with  time.  Computa- 
tions have  been  carried  out  showing  the  effect  of  damping  und  wn/w,  on  the  delay  r and  amplitude 
function  ft  of  the  response  curve,  it  is  found  that  for  larger  values  of  s , r reaches  a constant  value  less 
than  those  obtained  for  smaller  Also,  the  time  it  takes  to  reach  this  constant  value  is  much 
shorter  for  huger  values  of  £.  For  fixed  oJn/ojf.  decreasing  f increases  6,  while  for  fixed  f,  increasing 
u>n/tof  increases  6.  Results  for  input  loads  expressed  as  ramp  and  triangular  functions  are  also  pre- 
sented. For  the  triangular  input,  the  results  are  quite  similar  to  those  for  the  sinusoidal  case;  the  most 
significant  difference  being  the  smaller  values  obtained  in  the  response  curves. 
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FIG.  13:  RESPONSE  TO  A STEP  MODULATED  INPUT 
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DECAY  TO  A PULSE  MODULATED  INPUT  FOR  o>  /wf  = 0.5,  $ = 0.02 


DECAY  TO  A PULSE  MODULATED  INPUT 


DECAY  TO  A PULSE  MODULATED  INPUT 
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FIG.  20:  DECAY  TO  A PULSE  MODULATED  INPUT  FOR  c ojcot  = 1,  t = 0.04 


DECAY  TO  A PULSE  MODULATED  INPUT 


FiG.  22:  DECAY  TO  A PULSE  MODULATED  INPUT  FOR  w„/wf  = 0,  f = 0.08 


DECAY  TO  A PULSE  MODULATED  INPUT  FOR  a>Ja><  = 0.5, 


FIG.  24:  DECAY  TO  A PULSE  MODULATED  INPUT 


ODULATED I 


FIG.  26:  MEAN  SQUARE  RESPONSE 
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FfG.  23:  MEAN  SQUARE  RESPONSE  FOR  CASE  I:  t = 0.02, 


FIG.  30:  MEAN  SQUARE  RESPONSE  FOR  CASE  I : £ = 0.02,  VB  = 40,  to, 
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MEAN  SQUARE  RESPONSE  FOR  CASEI:  f = 0.02, 


FIG.  35:  MEAN  SQUARE  RESPONSE  FOR  CASE  I:  f = 0.02, 
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FIG.  38;  VARIATION  OF  AMPLITUDE  FUNCTION  5 WITH  DURATION  OF  APPLIED  FORCE 
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VARIATION  OF  AMPLITUDE  FUNCTION  5 WITH  DURATION  OF  APPLIED  FORCE 


FIG.  43:  VARIATION  OF  AMPLITUDE  FUNCTION  8 WITH  DURATION  OF  APPLIED  FORCE 


VARfATION  OF  AMPLITUDE  FUNCTION  5 WITH  DURATION  OF  APPLIED  FORCE 


VAR  I AT  JON  OF  AMPLITUDE  FUNCTION  5 WITH  DURATION  OF  APPLIED  FORCE 


EFFECT  OF  Lx'  ON  MEAN  SQUARE  RESPONSE  FOR  CASE  II 


MEAN  SQUARE  RESPONSE  FOR  CASE  II  AND  III: 


CASE  I!  AND  III:  ? = 0.02,  wn/w, 


MEAN  SQUARE  RESPONSE  FOR  CASE  II  AND  III:  f = 0.02,  ton/w 


FIG.  52:  MEAN  SQUARE  RESPONSE  FOR  CASE  II  AND  III:  f = 0.02 


FIG.  54;  MEAN  SQUARE  RESPONSE  FOR  CASE  If  AND  HI;  f = 0.02, 


ARE  RESPONSE  FOR  CASE  il  AND  III:  f = 0.02 


FIG.  57:  MEAN  SQUARE  RESPONSE  FOR  CASE  II  AND  III:  t = 0.02,  cjJoj. 


